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We consider the undamped nonlinear Schrodinger equation driven by a periodic external force. 
Classes of travelling solitons and multisoliton complexes are obtained by the numerical continuation 
in the parameter space. Two previously known stationary solitons and two newly found localised 
solutions are used as the starting points for the continuation. 

We show that there are two families of stable solitons: one family is stable for sufficiently low 
velocities while solitons from the second family stabilise when travelling faster than a certain critical 
speed. The stable solitons of the former family can also form stably travelling bound states. 
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I. INTRODUCTION 

The damped nonlinear Schrodinger equation driven by 
a time-periodic external force, 



lUt 



+ 2|upu + 5u = ae 



nt 



i^u, 



(la) 



and its parametrically driven counterpart model two 
fundamental energy supply mechanisms in a nearly- 
conservative spatially distributed system. While the un- 
perturbed Schrodinger is an archetypal equation for the 
slowly varying envelope of a group of dispersive waves, 
the damped-driven equations arise whenever the resonant 
forcing of small amplitude is used to compensate weak 
dissipative losses. 

The simplest (and perhaps the most visually appeal- 
ing) realisation of Eq. p^ is that of the amplitude equa- 
tion for a strongly coupled pendulum array with the hori- 
zontal sinusoidal driving JL], taken in its continuum limit. 
Here a and fi are the driving strength and driving fre- 
quency, respectively; 5 is the detuning of the driving fre- 
quency from the continuum of linear waves in the array, 
and j5 is the damping coefficient. 

The array of torsionally coupled pendula can serve as 
a prototype model for the whole variety of systems in 
condensed matter physics. Accordingly, Eq. ljlal) was em- 
ployed to study systems as diverse as the ac-driven long 
Josephson junctions and charge-density-wave conduc- 
tors with external electric field [3| ; double-layer quantum 
Hall (pseudo)ferromagnets B and easy- axis ferromagnets 
in a rotating magnetic field [5| . Eq JTal arises in the the- 
ory of rf-driven waves in plasma [a, |7[ and shear flows 
in nematic liquid crystals [SJ ; the same equation governs 
the amplitude of the slowly varying Tr-mode in the forced 
Fermi-Pasta-Ulam lattice [3] . 
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A closely related equation is the one with the spatially 
periodic forcing. 



iut + Uxx + 2|u| u + 5u = ae - i/3u, (lb) 

and, more generally, the one driven by the harmonic wave 
flMl: 



lUt 



■2\u\^u- 



- 5u = ae^(^^+"*) 



ijiu. (Ic) 



A discrete version of Ea. dlbp describes an array of 
coupled-waveguide resonators excited by a driving field 
[13J whereas Eq.([Tc]) models pulse propagation in an 
asymmetric twin-core optical fiber [10[ . 

Equation pc|) includes ([Ta)) and (jlbl) as particular 
cases. The transformation 

M(x,t) = *(A:,t)e*(^^+"*), X = x-2Kt 



takes Px)) to 



i^t + "^xx + 2|*p* - K^* = a - i/3^. 



(2) 



with K^ = K^ + Vt — 5. The equation in this form has a 
history of applications of its own — in particular, in the 
physics of optical cavities. Ori gina lly, it was introduced 
as the Lugiato - Levefer model [IJ of the diffractive cav- 
ity driven by a plane- wave stationary beam. Later it was 
employed to describe a synchronously pumped ring laser 
with a nonlinear dispersive fiber [1^, [l^ . More recently 
the same equation was shown to govern the envelopes of 
short baroclinic Rossby waves in the two-layer model of 
the atmosphere, or the ocean [17]. 

Equation ^ has undergone an extensive mathematical 
analysis. To pics covered included existence [IS. 1 2 01 l2ll|. 
stability [20, [221 and bifurcation [7, 1£| of nonpropagat- 



ing solitons and their bound states la, [23h25] ; statistical 
mechanics of soliton creation and annihilation 26]; soli- 
ton autoresonance phenomena [12, [23] ; regular [23 and 
chaotic [2^ attractors on finite spatial intervals. Here 



and below we use the word "soliton" simply as a syn- 
onyiTL for "localised travelling wave" . 

The recent paper [30| studied solitons of the undamped 
(/3 = 0) equation (^ travelling with constant or oscillat- 
ing velocities. Summarising results of their direct numer- 
ical simulations of Eq. ^ , the authors formulated an em- 
pirical stability criterion of the soliton against small and 
large perturbations. So far, this criterion has not been 
given any mathematical proof or physical justification. 
Despite being tested on a variety of initial conditions, it 
still has the status of conjecture. 

In order to verify the validity of the empirical stabil- 
ity criterion at least for infinitesimal perturbations, one 
needs to have the travelling soliton existence and lin- 
earised stability domains accurately demarcated. The 
classification of bifurcations occurring when stability is 
lost would also be a useful step towards the justification 
of the criterion. This is what we shall concern ourselves 
with in this paper. 

Here, we study travelling solitons of Eq.® by path- 
following them in the parameter space. One advantage 
of this approach over simulations is that it furnishes all 
soliton solutions moving with a given velocity — all sta- 
ble and all unstable. This, in turn, allows one to under- 
stand the actual mechanisms and details of the soliton 
transformations. 

The outline of this paper is as follows. In the next 
section, we give a brief classification of space- and time- 
independent solutions of Eq.© which may serve as the 
backgrounds for the solitons. In particular, we show that 
there is only one stable background and determine the 
value of the limit speed of the soliton propagating over 
it. In section Hill we describe insights one can draw from 
the analysis of the eigenvalues of the symplectic linearised 
operator and its hermitian counterpart. These pertain to 
the stability and bifurcation of the solitons. 

In section IIVI we present four nonpropagating directly 
driven solitons. Two of these have already been available 
in literature while the other two have not been known 
before. In sections |V] and IVIl we report on the contin- 
uation of these stationary solitons to nonzero velocities. 
Our results on the existence and stability of the trav- 
elling solitons and their complexes, are summarised in 
section IVlTl In particular, FiglH] gives a chart of "stable" 
velocities for each value of the driving strength. 



II. FLAT SOLUTIONS 



Assuming that k^ > and defining t' — K^t, x' = kX, 
and '5 = Ktp, equation ^ becomes 

#t' + ipx'x' + 2|?/'PV' -Tp ^ -h- ijip, 

where h = —a/n^, 7 — 13/ k^. (In what follows, we omit 
primes above x and t for notational convenience.) 

In this paper we study the above equation with zero 
damping: 7 = 0. Without loss of generality we can as- 
sume that ft, > 0. Since we shall be concerned with soli- 



tons travelling at nonzero velocities, it is convenient to 
transform the equation to a co-moving frame: 

#t-iyi/'? + V^?e + 2|V'lV-V' = -^, (3) 

where i_ = x — Vt. 

Flat solutions are roots of the cubic equation 

2\^\^~i:^~h- (4) 

these have been classified in [13]. If < ft < (2/27) ^/^ 
there are 3 roots, of which two (f/'i and ^2) are positive, 
and one (i/'s) is negative. Here f/;! < ^ < V'i < ^ < V'i < 
|. If ft > (2/27)^'^, there is only one (negative) solution 
V-s, with Vi > |. 

Let V'o denote a root of equation @ — one of the three 
roots ?Ai, -02 and f/'a- The value V'o does not depend on 
V: the flat solution has the same form in any frame of 
reference. However the spectrum of small perturbations 
of the flat solution does include a dependence on V. Let- 
ting tp = tpQ -\- [u{^) + iv{£))\e^^ in dS]), linearising in u 
and u, and, finally, taking u,v oc e^^^ , we obtain 



{X-ikVf 



-{k' + a'){e+b^), 



where we have introduced 



a = yr^ejg, b = ^1-2V;2. 



(5) 



(6) 



To determine whether i/jq can serve as a background 
to a stationary localised solution of ([3]) , consider a time- 
independent perturbation — that is, set A = 0: 



k^V^ = {k^+a^){k^+b''). 



(7) 



The only flat solution that is a priori unsuitable as a 
background for localised solutions is such -0o whose as- 
sociated quadratic equation ([7]) has two nonnegative real 
roots, (fc2)i > and {k^)2 > 0. 

It is not difficult to check that the negative solution 
'03 has two nonnegative roots for any choice of h and V. 
This disqualifies ips as a possible soliton background. We 
also conclude that travelling solitons may not exist for h 
greater than (2/27)^/2^ 

Next, ii V < c, where 

c = a + b, (8) 

the smaller positive solution ipi will have either two com- 
plex or two negative roots (fc^)i,2, whereas for velocities 
greater than c, both roots are nonnegative. Hence the 
"01 solution can serve as a background only for V < c. 
When V < b— a, the decay to the background is mono- 
tonic (both roots are negative), while when V > b — a, 
the decay is by ondulation (the roots are complex) . This 
flat solution admits a simple explicit expression: 



2 fa 27r 



where 



arccos — 




— < a < TT. 
2 - - 



Finally, the larger positive solution ■02 has two real 
roots of opposite signs (for aU V and < h < (2/27)^/2)^ 
This flat solution may also serve as a soliton background. 

Next, one can readily check that a flat solution ipo is 
stable if V'o ^ h- Therefore, even if there are solitons 
asymptotic to the flat solution ■!/;2 as x — >■ cx) or x — ^ — oo, 
these will be of little physical interest as the background 
■02 is always unstable. 

In summary, only the small positive flat solution (the 
one with ipQ < g) is stable. It may serve as a background 
for solitons only iiV<c; that is, the soliton propagation 
speed is limited by c. 

The inequality V < c limiting the soliton propagation 
speed, has a simple physical interpretation. Indeed, one 
can easily check that c gives the lower bound for the phase 
velocity of radiation waves [in the original (x, t) reference 
frame]. Therefore, a soliton travelling faster than c would 
be exciting resonant radiation. This is inconsistent with 
the asymptotic behaviour -^a, — >■ as |a;| — >■ cxd; neither 
could it be reconciled with the energy conservation. 



III. INSIGHTS FROM LINEARISATION 

Travelling wave solutions depend on x and t only in 
combination ^ = x~Vt. For these, the partial differential 
equation ^ reduces to an ordinary differential equation 



iVip^ + ip^^ + 2|V'PV' -ip^ ~h. 



(9) 



It is this equation that we will be solving numerically in 
the following sections. 

Let ipsiO be a localised solution of ([9]). In order to 
represent results of continuation graphically, we will need 
to characterise the function i/jsiO by a single value. A 
convenient choice for such a bifurcation measure is the 
momentum integral 



P = 



{iplip - ip^^*)d^. 



(10) 



One advantage of this choice is that the momentum is an 
integral of motion for equation ([3]); hence P is a physi- 
cally meaningful characteristic of solutions. Another use- 
ful property of the momentum is that in some cases its 
extrema mark the change of the soliton stability proper- 
ties (see below). 



A. The hermitian and symplectic operator 

Many aspects of the soliton's bifurcation diagram can 
be explained simply by the behaviour of the eigenvalues 
of the operator of linearisation about the travelling- wave 



solution in question. Therefore, before proceeding to the 
numerical continuation of travelling waves, we introduce 
the linearised operator and discuss some of its properties. 
Consider a perturbation of the solution of Eq.® of 
the form ip — tpg + [u{£,) + iv{£,)]e^'^ , with small u and v. 
Substituting -0 in Eq. ^ and linearising in u and v, we 
get a symplectic eigenvalue problem 



Hy = XJy. 
Here y is a two-component vector- function 



(11) 



m 



and H is a hermitian differential operator acting on such 
functions: 



n 



-dl 



h 1 - 2(37^2 + 12) ^vd^ - 47^I 

Vd^ - Am -dl + 1 - 2(3X2 + 7^2) 



with TZ and T denoting the real and imaginary part of 
the solution V's (C) ■ ips — Ti- + H- Finally, J is a constant 
skew-symmetric matrix 



J = 



-1 

1 



Assume that ■0s (C) is a localised solution decaying to 
"00 as X — >■ ±cxD, where 0q < g. The continuous spec- 
trum of the hermitian operator % occupies the posi- 
tive real axis with a gap separating it from the origin: 
E > Eq > Q. Discrete eigenvalues En satisfy P„ < Pq. 
On the other hand, the continuous spectrum of the sym- 
plectic eigenvalues (that is, the continuous spectrum of 
the operator J~^'H) occupies the imaginary axis of A out- 
side the gap (—iujQ, iioo). The gap width here is given by 



Wo 



(H 



a^){ki+b^)-Vko>0, 



(12) 



where fco is the positive root of the bicubic equation 

Discrete eigenvalues of the operator J^^H may include 
pairs of opposite real values A = ±p; pure imaginary 
pairs A = zbiw, with < cj < loq; and, finally, complex 
quadruplets A = ±p ± icj. 

We routinely evaluate the spectrum of symplectic 
eigenvalues as we continue localised solutions in V. If 
there is at least one eigenvalue A with ReA > 0, the solu- 
tion 0s is considered linearly unstable. Otherwise (that 
is, if all eigenvalues have ReA < 0), the solution is deemed 
linearly stable. 

B. Zero eigenvalues 

While the eigenvalues of the operator J~^'H (that 
is, the eigenvalues of the symplectic eigenvalue problem 



(fTTjl ) determine stability or instability of the solution tpg, 
the eigenvalues of the operator H are significant for the 
continuability of this solution. Of particular importance 
are its zero eigenvalues. 

At a generic point V, the operator H has only one 
zero eigenvalue, with the translational eigenvector ^Pj = 
{TZ^,I^). This is due to the fact that the stationary equa- 
tion ^ has only one continuous symmetry. For a given 
V, the solution ipsiS,) is a member of a one-parameter 
family of solutions ipsiS, — Q), where 9 is an arbitrary 
translation. 

On the other hand, the nonhermitian operator J^^H 
has two zero eigenvalues at a generic point. The reason 
is that the equation ([3]) as well as its linearisation, are 
hamiltonian systems. Real and imaginary eigenvalues of 
operators which generate hamiltonian flows always come 
in pairs: If /i is an eigenvalue, so is — /i [33|. The two 
zero eigenvalues of the operator J~^H reflect the fact 
that the function ips{£,)^ considered as a solution of the 
partial differential equation ([3]), is a member of a two- 
parameter family. One parameter is the translation; the 
other one is the velocity V. 

For generic V, the repeated zero eigenvalue of J^^H 
is defective: there is only one eigenvector \I^j associated 
with it. There is also a generalised eigenvector Vty, where 



dV dVj ■ 

This vector-function is not an eigenvector of J^^TL; in- 
stead, differentiating © in F one checks that z — 'i'v 
satisfies the nonhomogeneous equation 



nz = -.m. 



(13) 



[That is, ^v is an eigenvector of the square of the sym- 

plectic operator: {J^^'H)'^^v = 0.] 

As we continue in V, a pair of opposite pure-imaginary 
symplectic eigenvalues may collide at the origin on the 
A-plane and cross to the positive and negative real axis, 
respectively. The algebraic multiplicity of the eigenvalue 
A = increases from 2 to 4 at the point V — Vc, however 
if the hermitian operator H does not acquire the second 
eigenvalue E = at this point, the geometric multiplicity 
remains equal to 1 . The change of stability of the soliton 
solution does not affect its continuability, i.e. the soliton 
exists on either side oi V = Vc- In this case we have 
dP/dV = at the point where the stability changes [33 |. 

The continuation may be obstructed only when an- 
other (the second) eigenvalue of the operator H crosses 
through zero at y = Vc'- H^ — 0. If the corresponding 
eigenvector $ is not orthogonal to the vector-function 
J^ J in the right-hand side of equation p^ , its solution 
z — '^v will not be bounded. This implies a saddle- node 
bifurcation; the soliton solution ips cannot be continued 
beyond V = Vc- Note that although $ is an eigenvector 
of the symplectic operator J~^'H, the algebraic multiplic- 
ity of the symplectic eigenvalue remains equal to 2 in this 
case. 



Assume now that the eigenvector $ is orthogonal to 
J^j. This may happen if the soliton solution ipg of equa- 
tion ^ with y = Vc is a member of a iwo-parameter 
family of solutions "ijjs = 4's{^ — ^;x)j with x equal to 
some xo- Here we assume that each member of the fam- 
ily 4's{£, — (^\x) is a solution of Eq.® — with the same 



V ^Vc- Then $ is given by ^> ^ 
a root of the equation 



d-^/dx 



X=Xo 



Fix) = 0, 



where 



F{x)^ {^^,.m^)dt 



If Xo is 



(14a) 



(14b) 



the vectors ^^ and J^j will be orthogonal which, in turn, 
will imply that a bounded solution ^y of the equation 
(fT3)) exists. [In Ea. (ll4bp ( , ) stands for the R^ scalar 
product: (a, 6) = aibi + 0262.] In this case the value Vb 
is not a turning point; the soliton solution ip^ exists on 
both sides oi V = Vq. The algebraic multiplicity of the 
zero symplectic eigenvalue increases at the point V = Ve- 
in fact from the hamiltonian property it follows that it 
increases up to 4 (rather than 3). 

Recalling the definition of the momentum integral (jlOp 
and writing it in terms of the real and imaginary part of 
Tps, equation (J14p becomes simply 



dP 

dx 



= 0. 



X=Xa 



This condition ensures that a two-parameter family of 
solutions Vs(a; — 0;x), existing at the velocity V = Vq, 
has a one-parameter subfamily tps{x — 6;xo) continuable 

to y ^ K m. 



IV. NON-PROPAGATING SOLITONS 

A. Simple solitons 

The ordinary differential equation Q with V = 0, 

V-., + 2|7/-| V - ^ = -/i, (15) 

has two real- valued localised solutions, ^+ and ^_. 
These are given by explicit formulas J22| : 



■0±(a;) ^ipo 



2 sinh^ 



1 ± cosh /3 cosh( Ax) 



(16) 



where the parameter /3 (0 < /3 < 00) is in one-to-one 
correspondence with the driving strength h: 

V2cosh^/3 



(l + 2cosh^;3)3/2 
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FIG. 1. Stationary ?/)+ and ■;/)_ solitons 



As h increases from to ■\/2/27 « 0.2722, /3 decreases 
from infinity to zero. (Hence < /i < 0.2722 is the 
domain of existence of the two sohtons.) The asymptotic 
value "00 E^nd inverse width A are also expressible through 

/3: 



-00 



1 



\/2 a/T+TcosI?^' 



A 



y2sinh^ 
VT+Tcosh^' 



(Note that the asymptotic value ■00 corresponds to the 
stable background, denoted Vi in SecHIl) 

The stationary soliton ■)/)+ has a positive eigenvalue in 
the spectrum of the linearised operator (jlll) : hence the 
-0+ is unstable for all h for which it exists [22]. The 
spectrum of the stationary soliton 0_ with small h in- 
cludes two discrete eigenvalues Ai,2 — Jwi,2, ^\,2 > — 
and their negative-imaginary counterparts. As h grows 
to 0.07749, Ai and \i approach each other, collide and 
acquire real parts of the opposite sign. This is a hamil- 
tonian Hopf bifurcation. For h > 0.07749, the soliton -0- 
is prone to the oscillatory instability |22l |. 

When a damping term is added to the equation, the 
two stationary solitons ip+ and ip- persist and can form 
a variety of multisoliton bound states, or complexes 
[m [23l - [25{ . In the next subsection, we show that un- 
damped directly driven solitons can also form stationary 
complexes. Some of these complexes are bound so tightly 
that the solution represents a single entity. To distin- 
guish these objects from the solitons -0+ and ?A_ , we will 
be referring to the ^+ and ■0- as the simple solitons. 



B. The twist solitons 

In addition to the two simple solitons expressible in 
elementary functions, the stationary equation (|15p has 
two localised solutions that cannot be constructed an- 
alytically. Unless h is extremely small, each of these 



two solutions has the form of a single entity [Figl3Ka,b)] 
— a soliton whose phase does not stay constant but 
grows, monotonically, as x changes from large negative 
to large positive values. When visualised in the three- 
dimensional (a;. Re?/;, Im0)-space, it looks like a twisted 
ribbon (twisted by 360°); hence we will be calling these 
two solutions simply "twists". For the reason that will 
become obvious in the paragraph following the next one, 
we denote the two solutions "0X2 and ipxa , respectively. 

The twist solitons were previously encountered in the 
parametrically driven (undamped) nonlinear Schrodinger 
equation [32| . For each h, the parametrically driven twist 
is a member of a two-parameter family of stationary two- 
soliton solutions. The first parameter is the overall trans- 
lation of the complex; the second one is the separation 
distance between the two bound solitons. The twist cor- 
responds to a very small separation, where the two sim- 
ple solitons bind to form a single entity. (The resulting 
object does not have even a slightest reminiscence of a 
two-soliton state; without knowing the whole family, the 
relation would hardly be possible to guess.) 

The two simple solitons, ^0+ and 0_, detach from the 
[/(l)-symmetric family of solitons of the unperturbed 
nonlinear Schrodinger at /i = [21[. The two twist so- 
lutions of (|15l) also hail from the solitons of the unper- 
turbed equation; however this time the relation is more 
complicated. Reducing /i, the two solutions transform 
into complexes of well-separated solitons [Figl2Jc,d)]. 
Namely, one of the two twist solutions becomes a complex 
of two solitons: 

ipT2 -^ e^'''/'^sech{x + xo) + e-^'''^'^sech{x - xq), 

where xq -^ oo as h -^ 0. The other twist continues to a 
complex of three unperturbed solitons: 



0T3 — >■ i sech(a; + xo) — sechx — i sech(a 



xo) 



and again, the separation xq grows without bound as 
/i — >■ 0. The "full names" of the two twists, 0t2 and ipT3, 
were coined to reflect this multisoliton ancestry. 

Despite being quiescent, nonpropagating objects, the 
twists carry nonzero momentum. Since equation (jlSp 
is invariant under the space inversion, the twist soliton 
with momentum P has a partner with momentum — P 
which is obtained by changing x -> — a;. This trans- 
formation leaves the absolute value of ip{x) intact but 
changes the sign of the phase derivative, {d/dx)a,Tgil){x). 
By analogy with the right-hand rule of circular motion, 
the twist whose phase decreases as x grows from — oo 
to -|-oo [that is, the trajectory on the (Re0,lm'0) phase 
plane is traced clockwise], will be called right-handed. 
The twist with the increasing phase (i.e. with a trajec- 
tory traced counter-clockwise) will be called left-handed. 
One can readily verify that the left-handed twist has a 
positive momentum, whereas the right-handedness im- 
plies P < 0. 
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FIG. 2. (a,b): The two nonpropagating twist solutions for h ^ 1. (Here h = 0.2). (c,d): The corresponding quiescent solutions 
when continued to an exponentially small h. (Here h — 2.15 x 10~^). All twist solutions shown in these figures are left-handed. 



Consider some particular value of the driving strength, 
h — ho- Unlike the twist solution in the parametrically 
driven NLS, the directly driven twist with h = hi^ is a, 
member of a one-parameter (rather than two-parameter) 
family of solutions. (The only free parameter is the trans- 
lation, — oo < < oo, whereas the intersoliton separation 
X is fixed by /i.) This can be concluded from the fact 
that the corresponding operator % has only one, trans- 
lational, zero eigenvalue. Had the twist been a member 
of a family of solutions parametrised by two continuous 
parameters, say 6 and x, the operator % would have had 
an additional zero eigenvalue with the eigenvector ^i^. 

Letting ij) — xi-\- 1x2, the stationary equation (jlSp can 
be written as a classical mechanical system on the plane, 
with the Lagrangian 



L^-(il+xl) 



1 



{xl 



-Af 



1 



{x\ 



hxi. 



The existence of a one-parameter family of homoclinic 
orbits X = x^(t), where x = (a;i,a;2), would imply that 
the above system has the second integral of motion, in 
addition to the energy. However, equation (jlSp is known 
not to have any additional conserved quantities |34| . 

Finally, we need to comment on the stability of the two 
twist solutions. When h is equal to and the two solu- 
tions represent a doublet and a triplet of infinitely sepa- 
rated solitons of the unperturbed nonlinear Schrodinger, 
the symplectic spectrum includes 8 and 12 zero eigen- 
values, respectively. When h is small nonzero, only two 
eigenvalues remain at the origin in each case. In addi- 
tion, the spectrum of the '0t2 twist includes a complex 
quadruplet ±A, ±A* and a pair of opposite pure imagi- 
nary eigenvalues. As h is increased, the imaginary pair 
collides with another imaginary pair emerging from the 
continuum, producing the second complex quadruplet. 
The spectrum of the ^t3 twist includes two complex 
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FIG. 3. As 1^ ^^ c, the tp+ solitons (for all h) and tp- soli- 
tons (for small h) approach linear waves with slowly decaying 
envelopes. Shown is the i/)+ solution with V close to c. (In 
this plot, h = 0.01; the corresponding c = 1.9996.) The ip- 
solutions with V close to c have a similar shape. 



■0+ . For a sequence of h sampling the interval (0, •y/2/27) , 
the branch starting at ■0+ was path followed all the way 
to V = c, where c is given by Eq.([8l). As V increases, 
the amplitude of the solution decreases while the width 
grows. A typical solution with V close to c is shown in 
Fig El As y — > c, the momentum P tends to zero. 

The resulting P{V) diagram is shown in FigHJ^a). For 
each h, the unstable stationary 0+ soliton remains un- 
stable when travelling sufficiently slow. The instability is 
due to a real eigenvalue A > of the linearised operator 

As V grows, the unstable eigenvalue moves towards the 
origin along the real axis. Eventually, as the momentum 
P reaches its maximum, the positive eigenvalue A collides 
with its opposite partner A' = —A, after which both real 
eigenvalues move onto the imaginary axis and the soliton 
acquires stability. The soliton remains stable all the way 
from the point Vc, where the momentum is maximum, to 
the value V — c where P = and the soliton ceases to 
exist. 

The resulting P(V) dependence shows a remarkable 
similarity to the P(V) diagram |32l] for the parametrically 
driven nonlinear Schrodinger, 



quadruplets and a pair of pure imaginary eigenvalues; 
this arrangement remains in place for all h, from very 
small to h = -^2/27. The bottom line is that both twist 
solutions are unstable for all h; the instability is always 
of the oscillatory type. 



V. NUMERICAL CONTINUATION OF SIMPLE 
SOLITONS 

A. The travelling 1/1+ soliton 

Travelling solitons are sought as solutions of the ordi- 
nary differential equation ^ under the boundary condi- 
tions Vc — ?► as 1^1 — >■ 00. 

We begin with the continuation of the quiescent soliton 



iijjt - iVil)^ + 044 + 2|V'PV' - V' = hil)* 



(17) 



The "parametrically driven" diagram is reproduced in 
FiglUb) for the sake of comparison. One should keep in 
mind here that the notation used for the parametrically 
driven solitons is opposite to the notation employed in 
the externally driven situation. Thus, the parametrically 
driven stationary {V = Q) soliton with a positive sym- 
plectic eigenvalue in its spectrum is denoted V_ (and 
not Tp+ as its externally driven counterpart). On the 
other hand, the parametrically driven stationary soliton 
denoted 0+ is stable for sufficiently small h (like the ex- 
ternally driven soliton '(/'-)• For this reason, the objects 
featuring P{V) diagrams similar to those of our exter- 
nally driven solitons ip^^ are the parametrically driven 
solitons ip-. 



B. The travelling t/j. soliton; h < 0.06 

In the case of the 0_ solitons, there are two charac- 
teristic scenarios. When h lies between and 0.06, the 
soliton ip- exists for all V between and c. As V is in- 
creased from zero, the momentum P grows from P = Q 
and reaches its maximum at some point Vc, < Vc < c. 
As V is changed from V to c, the momentum decays to 
zero [see Figllja)]. On the other hand, when h equals 
0.06 or lies above this value, the curve P{V) does not 
exhibit a point of maximum. 

Consider, first, the case h < 0.06. The transformation 
scenario here is similar to the case of the soliton ip+ ; see 
FiglD What makes the bifurcation curves for the 0+ 



r 



and ip^ solitons different, is the stability properties of 
the two solutions. Unlike the '0+ solution, the stationary 
■0_ soliton with h < 0.07749 is stable and its stability 
persists when it is continued to small nonzero velocities. 
As V grows to the value Vc where the momentum reaches 
its maximum, two opposite pure imaginary eigenvalues 
collide at the origin on the (ReA, ImA) plane and cross 
to the positive and negative real axis, respectively. For 
the driving strengths h < 0.055, this implies the loss of 
stability. 

As for the interval 0.0551 < h < 0.06, here the instabil- 
ity sets in earlier, as V reaches some V = Vq (where Vo < 
Vc). At the point V = Vo, two pairs of pure imaginary 
eigenvalues collide and produce a quadruplet of complex 
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FIG. 4. (a) The momentum of the 1/)+, ■)/'_ sohtons continued to positive velocities. Decimal fractions attached to branches 
label the corresponding values of /i, with the superscripts + and — indicating the ■)/'+ and i/i- solitons. (For example, 0.1"^ 
marks the branch emanating from the stationary i])^ soliton with /i = 0.1.) The only two branches that have not been labelled 
are the i/>+ and ?/)_ branches with h — 0.01; the curve just above ft = is 0.01~ and the curve just below the h — branch is 
O.Ol""". Solid curves mark stable and dashed ones unstable branches, (b) The corresponding P{V) diagram for the travelling 
parametrically driven solitons from 32J. 



eigenvalues ±A,±A*. (Here A has a small real and fi- 
nite imaginary part.) This is a point of the hamiltonian 
Ho pf b ifurcation, associated with the oscillatory instabil- 
ity ^[l!]. As F is increased to Vi (where Vb < Fi < V^), 
two pairs of complex-conjugate A converge on the real 
axis, becoming two positive (Ai = A2 > 0) and two neg- 
ative (— Ai = — A2) eigenvalues. Finally, when V crosses 
through Vc, the eigenvalues Ai and — Ai move on to the 
imaginary axis. The soliton does not restabilise at this 
point though; the real pair ±A2 persists for all V > Vc- 

The bifurcation values Vq and Vi are, naturally, func- 
tions of h. The value Vb decreases (and Vi increases) as 
h is increased from 0.0551. Eventually, when h reaches 
0.07749, Vo reaches zero. It is interesting to note that 
there is a gap between Vi and Vc for all h. Therefore 
the oscillatory and nonoscillatory instability coexist for 
no V; for smaller V {Vq < V < Vi) the instability is 
oscillatory whereas for larger V {V > Vi) the instability 
has a monotonic growth. 

Finally, it is appropriate to mention here that the bi- 
furcation curve for the ip^ solitons with small h < 0.06 
has the same form as the P{V) dependence for the small- 
h parametrically driven solitons (more specifically, para- 
metrically driven ip-plus solitons) — see FigUJ^b). 



C. The travelling V- soliton; h > 0.06 

The P{V) graphs for h > 0.06 are qualitatively differ- 
ent from the small- /i bifurcation curves. For these larger 
h, the bifurcation curve emanating from the origin on 



the {V, P)-plane turns back at some V = Vmax, with 
the derivative dP/dV remaining strictly positive for all 

V < Knax. 

For h in the interval 0.06 < h < 0.25, the P{V) curve 
crosses the P-axis [FigHJa), FigjS]. The solution arising 
at the point 1/ = is nothing but the '0T2 twist soliton, 
shown in Fig 121 a). 

As we continue this branch to the V < 0-region, the 
twist transforms into a complex of two well-separated ^p- 
solitons. The P{V) curve makes one more turn and even- 
tually returns to the origin on the (y,P)-plane (FiglSj. 
As V and P approach zero, the distance between the 
solitons in the complex tends to infinity. 

An interesting scenario arises when h is greater or 
equal than 0.25. Here, as V grows from zero, the soli- 
ton ?/;_ gradually transforms into a three-soliton complex 

-0(^ i_). The branch turns back towards ^ = but does 

not cross the P-axis. Instead of continuing to negative V, 
the branch reapproaches the origin in the (V, P) plane, 
remaining in the positive (V, P) quadrant at all times. 
The ingoing path is almost coincident with the outgoing 
trajectory; as a result, the branch forms a lasso- looking 
loop [Figgi;a)]. 

Turning to the stability properties of solutions along 
the branch continued from '0-, we start with a short in- 
terval 0.06 <h< 0.07749. The movements of the stabil- 
ity eigenvalues along the section of the curve emanating 
from the origin on the {V, P) plane, are similar to the in- 
terval 0.055 < h < 0.06 that we discussed in the previous 
paragraph. The stationary ?/;_ soliton is stable and sta- 
bility persists for small V. As V reaches a certain Vq > 0, 




FIG. 5. The full P{V) bifurcation diagram for the tp- soli- 
ton with h = 0.2. Also shown is the continuation of the T3 
twist and the ip(++) branch. More solution branches can be 
obtained by the reflection V — >■ —V, P — >■ —P. All branches 
shown in this figure correspond to unstable solutions. 



a quadruplet of complex eigenvalues is born and oscil- 
latory instability sets in. Subsequently two pairs of the 
complex eigenvalues converge on the real axis, dissociate, 
then recombine and diverge to the complex plane again; 
a pair of opposite pure imaginary eigenvalues moves to 
the real axis and back — however, despite all this activity 
on the complex plane, the soliton solution never regains 
its stability. 

For larger h, h > 0.07749, the stationary ip^ soliton 
is unstable, with a complex quadruplet in its spectrum. 
As we continue in V, two pairs of opposite pure imagi- 
nary eigenvalues move on to the real axis, one after an- 
other. For h > 0.25, the resulting arrangement (two pairs 
of opposite real eigenvalues and a complex quadruplet) 
persists until the branch reaches the origin on the {V, P) 
plane. On the other hand, when h lies in the interval 
0.07749 < h < 0.25, the four real eigenvalues collide, 
pairwise, producing the second complex quadruplet at 
some point on the curve before it crosses the P axis in 
FigsHJa) and[5l Two complex quadruplets persist in the 
spectrum as we continue the curve further. Thus the un- 
stable stationary soliton V'- with h > 0.07749, remains 
unstable for all V. 



VI. NUMERICAL CONTINUATION OF THE 
TWIST SOLITON 

When 0.06 < h < 0.25, the branch resulting from the 
continuation of the stationary ijj- soliton turns back and 
crosses the P-axis; the point of crossing corresponds to 



the T2 twist solution. On the other hand, when h lies 
outside the (0.06,0.25) interval, the stationary T2 twist 
is disconnected from the stationary ■)/)_ soliton and can be 
used as a starting point for a new, independent, branch. 
Another new branch is seeded by the T3 solution. 

These additional branches of travelling solitons are 
traced in this section. 



A. Travelling twist T2 (h < 0.06) 

We start with the situation of small h: h < 0.06, and 
consider the T2 solution first. 

When the stationary ^t2 twist is path followed to pos- 
itive V, it transforms into a '4'{++) complex. At some 
point, the P{V) curve makes a U-turn [Fig. [S]^a)] and 
connects to the origin on the (V, P) plane. The entire 
positive- T^ branch is unstable. The stationary twist has 
a complex quadruplet in its spectrum; as the curve is 
continued beyond the turning point, the complex eigen- 
values converge, pairwise, on the positive and negative 
real axis. In addition, a pair of opposite pure imaginary 
eigenvalues moves onto the real axis as V passes through 
the point of maximum of the momentum in FigjSfa). 

As the curve approaches the origin, the distance be- 
tween the two solitons in the complex increases and be- 
comes infinite when V = P = 0. The spectrum becomes 
the spectrum of two infinitely separated '0-1- solitons, i.e. 
it includes two positive eigenvalues Ai « A2; their nega- 
tive counterparts — Ai ~ — A2; and four eigenvalues near 
the origin. 

Continuing the T2 twist in the negative- T^ direction, it 
transforms into a complex of two tp^ solitons. At some 
point along the curve, a quadruplet of complex eigen- 
values converges on the imaginary axis and the complex 
stabilises. (For the value h — 0.05 which was used to 
produce FigEfa), the stabilisation occurs at the point 

V = —0.45.) Continuing to larger negative V, the 
branch turns back; shortly after that (at V = —0.503 
for h — 0.05) the momentum reaches its minimum. Two 
opposite imaginary eigenvalues collide at this point and 
move onto the real axis; the solution loses its stability. 

When continued beyond the turning point and the 
point of minimum of momentum, the curve connects to 
the origin on the (V, P) plane (Fig. EKa)). As V, P -> 0, 
the distance between the two ip- solitons grows without 
bound. The two opposite real eigenvalues decay in ab- 
solute value but remain in the spectrum all the way to 

V = {). 

It is interesting to note a similarity between the bi- 
furcation diagram resulting from the continuation of the 
small-/i T2 twist in the externally driven NLS [FigjS] 
(a)] and the corresponding diagram in the parametrically 
driven case. The latter is reproduced, for convenience of 
comparison, in Fig|S] (b). In both cases the continuation 
of the twist solution to negative velocities gives rise to a 
stable complex of two stable solitons. 
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FIG. 6. (a) Continuation of the T2 and TS twist solutions in the case of small h. Also shown is the two-soliton branch which 
connects the origin to itself without intersection the vertical axis, (b) The continuation of the twist soliton in the case of the 
parametrically driven NLS equation (adapted from [32|). More solution branches can be obtained by the reflection V -^ —V ^ 
P ^ -P both in (a) and (b). 



B. Travelling twist T3, h <0.25 

FigsUJa) and \E\ also show the continuation of the T3 
twist soliton. The bifurcation diagrams obtained for h < 
0.06 and 0.06 < h < 0.25 are qualitatively similar. 

Continuing the stationary T3 to positive velocities, the 
solution transforms into a '!/'(h h) complex. If we, in- 
stead, continue to negative velocities, the twist trans- 
forms into a triplet of tp- solitons. Both V > and 
T^ < parts of the curve turn and connect to the origin 
on the (y, P) plane. As V and P approach the origin on 
either side, the distance between the three solitons bound 
in the complex grows without limit. 

The stationary T3 has two complex quadruplets in its 
spectrum; depending on h, both or one of these converge 
on the real axis as we continue it to F > and I^ < 0. 
Two opposite eigenvalues cross through A = at the 
extrema of P{V). Finally, as V and P approach the 
origin, the spectrum transforms into the union of spectra 
of three separate solitons. 



C. Travelling twists T2 and T3, h > 0.25 

Another parameter region where the continuation of 
the tjj- does not cross the P-axis, is h > 0.25. The result 
of the continuation of the two twist solutions is shown 
in FiglTl^a). The continuation of T2 to the negative ve- 
locities proceeds according to scenario similar to h = 0.2 

I 



D. Other branches 



and h = 0.05: the twist transforms into a complex of two 
solitons ip-. At some negative V the curve turns back 
and connects to the origin on the {V, P) plane, with the 
distance between the two solitons bound in the complex 
increasing without bound. The eigenvalues evolve ac- 
cordingly: two complex quadruplets in the spectrum of 
the stationary T2 persist for all F < 0, supplemented by 
a pair of real eigenvalues which arrive from the imaginary 
axis at the point of minimum of P{V). As V, P — > 0, the 
discrete spectrum becomes the union of the eigenvalues 
of two simple solitons. 

The continuation of r2 to positive V produces a less 
expected outcome. Instead of turning clockwise and con- 
necting to the origin as in FiglSfa), the curve turns coun- 
terclockwise and crosses through the P-axis once again. 
The solution arising at the point T^ = is nothing but 
the twist T3. Two complex quadruplets in the spectrum 
of T2 persist as it is continued to T3. 

The subsequent continuation produces a hook-shaped 
curve similar to the curve described in the previous para- 
graph and leading to the origin on the {V, P)-plane. The 
corresponding solution is a complex of three tp- soli- 
tons, shown in FiglTl^b). The third complex quadruplet 
emerges at some V before the turning point, and a pair of 
opposite real eigenvalues arrives from the imaginary axis 
at the point of minimum of the momentum. As V, P — > 0, 
the distance between the solitons grows to infinity and 
the spectrum approaches the union of the eigenvalues of 
three separate solitons '0_ . 



I 

listed above {ip±, '0T2 or ^ts)- The simplest of these 



It is appropriate to note that there are branches which 
do not originate on any of the four stationary solutions 
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emerge from the origin on the (V, P) plane as bound 
states of simple solitons with large separation. One 
branch of this sort arises for h > 0.06 (Fig. [SJ. It emerges 
from the origin as the '0(++) a-nd returns as the '(/'(++++) 
complex. The entire branch is unstable. 

Next, unlike in the parametrically driven NLS, the 
same pair of externally driven travelling solitons may 
bind at various distances. In particular, when h is smaller 
than 0.06, there is more than one bound state of two ip+ 
solitons and more than one complex of two '0-minuses. 

Figinfa) shows a branch ipr ) that emerges from the 

origin in the first quadrant of the (V, P) plane, describes 
a loop and re-enters the origin — this time as a V'(+-i-) 
branch. Note that for small V and P, the re-entering 
'0(++) branch is indistinguishable from the other '/'(++) 
branch — the one that continues from the twist solution. 
(In a similar way, the V^— > —V, P— > —P refiection of 

the V'( ) branch overlaps with the small- V, P section of 

the ^pt ) branch arriving from the twist.) All solutions 

constituting this branch are unstable. 



VII. CONCLUDING REMARKS 

In this paper, we studied stationary and moving soli- 
tons of the externally driven nonlinear Schrodinger equa- 
tion. 



i->pt + i>xx + SIV'PV' - V" = -h. 



(18) 



Our continuation results are summarised in FiglHlJa) 
which shows ranges of stable velocities for each value of 
the driving strength h. 

The notation tp^ and ip- in this figure is used for the 
travelling waves obtained by the continuation of the sta- 



tionary ip^ and V'- solitons, respectively. The travelling 
soliton preserves some similarity with its stationary an- 
cestor; this justifies the use of the same notation. 

The uppermost curve in this figure is given by 1^ = c{h) 
where c is the maximum velocity of the soliton propaga- 
tion, Eq. ([HI) . This curve serves as the upper bound of the 
travelling V'+ soliton existence domain. The dotted curve 
demarcates the existence domain of the travelling ip^ soli- 
ton. For h between and 0.06 it coincides with the V ~ c; 
for 0.06 < h < 0.2722 it is given by F = Kiax(^) where 
Vmax is the position of the turning point in Fig|4l[a). 

The area shaded in blue (light grey) gives the stability 
region of the soliton -0+ and the area shaded by purple 
(dark grey) is the ip- stability domain. Note that the 
blue and purple regions partially overlap: for small h, 
there is a range of "stable" velocities accessible to solitons 
of both families. The light (yellow) strip inside the purple 
(dark grey) region represents the stability domain of the 
bound state of two V'- solitons. 

As we cross the right-hand "vertical" boundary of the 
purple (dark grey) region, the tp- soliton loses its stabi- 
lity to an oscillatory mode. If we had damping in the 
system, the onset of instability would correspond to the 
Hopf bifurcation giving rise to a time-periodic solution. 
In the absence of damping, the oscillatory instability pro- 
duces an oscillatory structure with long but finite lifetime 
[3l| . These solitons with oscillating amplitude and width, 
travelling with oscillatory velocities, were observed in 
|30l ] . These are expected to exist to the right of the purple 
(dark grey) region. 

Where possible, we tried to emphasise the similar- 
ity of the arising bifurcation diagrams with the corre- 
sponding diagrams for the parametrically driven nonlin- 
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FIG. 8. (a) The chart of the stable one-soliton solutions of the externally driven nonlinear Schrodinger equation pSjI . Here h 
varies from to -\/2/27 « 0.2722. (b) The corresponding attractor chart for the parametrically driven soliton (adapted from 



ear Schrodinger equation: 



(19) 



FigHfb) reproduces the soliton attractor chart for 
Eq. P^ |32| . The structure of the stability regions in the 
two figures is remarkably similar. The slowly moving soli- 
tons in the purple- (dark grey-) tinted region inherit their 



stabiUty from the stationary solitons of the family which 
is stable for small h (the "ip- family in the externally- 
driven and the ■0-1- family in the parametrically-driven 
case). On the other hand, the solitons in the blue- (light 
grey-) shaded area are transonic (i.e. move close to c, 
velocity of the sound waves). Their stability is due to 
the proximity of the nonlinear Schrodinger equation to 
the KdV in the transonic limit 35] . 
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